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Field of the Invention 

This invention relates to flow fields for uniform distribution of fluids or their 
active components or properties to and from a target area. The flow field may be 
embodied in a flow field device such as a flow field plate or bipolar plate used for 
10 distribution of reactants to, and removal of products from, opposite sides of a catalyzed 
membrane in an electrochemical cell such as a fuel cell. 



Background of the Invention 

A number of references depict flow fields having serpentine channels wherein 
15 sequential segments of the channels are parallel, including: U.S. Pats. Nos. 4,686,159; 
4,988,583; 5,108,849; 5,252,410; 5,683,828; 5,750,281; 5,773,160; 5,846,668; 
5,858,567; 5,858,569; 5,945,232; 6,071,635 and 6,099,984. 

A number of references depict flow fields having multiple interleaved 
serpentine channels wherein sequential segments of each channel are parallel, 
20 including: U.S. Pats. Nos. 5,683,828; 5,750,281; 5,773,160; 5,804,326; 5,840,438; 
5,858,567; 5,998,055; 6,071,635 and 6,093,502. 

A number of references depict interdigitated flow fields, including: U.S. Pats, 
Nos. 5,252,410; 5,641,586 and 6,207,312. In an interdigitated flow field, channels 
having an inlet but no outlet alternate with channels having an outlet but no inlet. 
25 In addition, the use of a metal screen as a flow field has been taught, e.g. in U.S. 

Pat. Nos. 4,855,193; 5,798,187; 6,037,072 and 6,207,310. 

U.S. Pat. No. 5,922,485 depicts flow fields having serpentine channels 
composed of concentric circular segments, as well as straight-line serpentine channels. 
U.S. Pat. No. 5,686,199 depicts a series-parallel arrangement composed 
30 essentially of parallel segments. 
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U.S. Pat. No. 6,048,634 depicts flow field patterns wherein pairs of adjacent 
channels carry flow in opposite directions, including spiral patterns and serpentine 
patterns wherein sequential segments of the channels are parallel. 

U.S. Pat. Nos. 4,631,239 and 4,853,301 describe serpentine flow fields wherein 
5 sequential segments of the channels are parallel, where the segments are skewed 
relative to the boundaries of the bipolar plate and/or relative to the flow field on the 
opposite face of a bipolar plate. 

U.S. Pat. No. 4,292,379 describes flow fields on either side of a bipolar plate 
wherein the depth and/or separation of parallel channels are varied so as to create an 
10 uneven distribution that matches the uneven distribution created by the opposing face of 
the plate. 

U.S. Pat. No. 4,324,844 concerns an electrochemical cell that includes cooling 
fluid flow passages having varying surface area and spacing. 

15 Summary of the Invention 

Briefly, the present invention provides a fluid distribution assembly comprising 
a flow field device embodying a flow field and a fluid transport layer disposed between 
the flow field device and a target area, where, for at least one finite non-zero flow rate 
and at least one use rate of an active component or property of the fluid in the fluid 

20 transport layer, lateral flux of the active component or property varies by no more than 
35% through at least 90% of all overland portions of said fluid transport layer. In one 
embodiment, the flow field device comprises a flow field comprising a serpentine 
channel, comprising non-parallel sequential major segments. In a further embodiment, 
the angles between successive major segments of the serpentine channel vary 

25 progressively. 

In another aspect, the present invention provides a flow field device embodying 
a flow field comprising at least one serpentine channel wherein at least two sequential 
major segments of the channel are non-parallel. 

In another aspect, the present invention provides a flow field device embodying 

30 a flow field comprising a channel, where the spacing between analogous parts of 

sequential major segments of channel decreases monotonically with distance from the 
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inlet, or where land areas separating the major segments decrease in size monotonically 
with distance from the inlet. 

What has not been described in the art, and is provided by the present invention, 
is a flow field designed to provide uniform lateral flux through the fluid transport layer, 
5 and in particular by use of a "zig-zag" serpentine or "progressive" flow channel. 

In this application: 

"flow field" refers to a pattern of one or more channels embodied in a 
component of a fluid distribution system, which system allows ingress and egress of 
fluids to and from a target area; 
10 "target area" refers to an area having significant extent in at least two 

dimensions which is served by a fluid distribution system, such as the electrochemically 
active electrode area of an electrochemical device; 

"active area" refers to the area of a flow field overlaying and serving the target 

area; 

15 "active component" refers to a component of a fluid to be used at or in 

conjunction with the target area, e.g. the oxygen present in air, the hydrogen present in a 
reformate gas mixture, and the like; 

"active property" refers to a property of a fluid to be used at or in conjunction 
with the target area, e.g. the thermal energy content of a coolant, the solvating capacity 

20 of a solvent, and the like; 

"serpentine" refers to a pattern, such as the pattern of a channel in a flow field, 
comprising sequentially connected major segments which alternate in orientation, such 
as in orientation of flow direction, and which meet at turning points or are connected by 
turning segments; 

25 "land" or "land area" refer to area between channels or portions of channels of a 

flow field; 

"major segment" refers to a segment of a pattern, such as the pattern of a 
channel in a flow field, having a geometrical orientation distinct from that of major 
segments directly preceding or following, which is connected to major segments 
30 directly preceding or following either at tuming points or through relatively short 
turning segments; 
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"flux" refers to the transport of a fluid, such as a gas or liquid, which can be 
expressed in units of kg/s/m^, or the transport of a component of a fluid, e.g. the oxygen 
present in air, through a given area, which can be expressed in units of kg/s/m^ or the 
transport of a property of a fluid, e.g. thermal energy which can be expressed in units of 
5 watts/m^ 

"flow rate" refers to the transport of a fluid, such as a gas or liquid, or the 
transport of a component of a fluid, e.g. the oxygen present in air, and can be expressed 
in units of mass per unit time (e.g., kg/s) or volume at standard conditions per unit time 
(e.g., standard cubic centimeters per minute (seem) or standard liters per minute (slm)); 
10 "fluid transport layer" means a layer allowing fluid transport, typically a layer of 

porous or otherwise fluid-permeable structural material, but also including a gap 
CI maintained mechanically; 

iji "overland portions" of a fluid transport layer used with a flow field are portions 

■ y of the fluid transport layer that pass over a land area within the active area of the flow 

2 15 field, which excludes portions of the fluid transport layer that are over a flow field 

channel or portions not passing over the active area of the flow field; 
t^; "lateral flux" of fluid through a layer, such as a fluid transport layer, means flux 

jO within the layer and generally within the plane of the layer, as distinguished from flux 

M into or out of the layer that may be flux orthogonal to the plane of the layer; 

2 20 "flow field device" refers to a component of a fluid distribution system which 

embodies a flow field, typically a component of a fluid distribution system in an 
electrochemical cell, which is typically either i) a flow field plate or ii) a fluid transport 
layer that is sufficiently sturdy to hold the pattern of a flow field stamped, molded or 
cut therein, but which is more typically a flow field plate; which is typically a bipolar 
25 plate, which may be made of porous or more typically non-porous material, and which 
is typically made of electrically conductive material. 

It is an advantage of the present invention to provide flow fields and flow field 
devices capable of highly uniform distribution of fluids or their active components over 
a target area, which may provide more efficient, stable and durable function in an 
30 electrochemical cell. 
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Brief Description of the Drawings 

Fig. 1 depicts a computational fluid dynamic (CFD) model of the comparative 
flow field of Example IC, having a single serpentine channel composed of parallel 
sequential channels, with calculated gas velocities in m/sec represented according to the 
5 indicated gray scale. 

Fig. 2 is a schematic depiction of a "zig-zag" serpentine flow field according to 
the present invention depicting parameters discussed in the text. 

Fig. 3 depicts a CFD model of the "zig-zag" serpentine flow field according to 
the present invention of Example 2, with calculated gas velocities represented 
10 according to the indicated gray scale, which is the same gray scale as Fig. 1 . 

Fig. 4 is a graph of inlet-to-outlet pressure drop vs. mass flow rate for the flow 
O fields of Example 2 and comparative Example 1 , 

Fig. 5 is a graph of gas velocity crossing a land feature vs. lateral position along 
the land feature for each of Examples 2-4 and comparative Example 1 . 
ipl 15 Fig. 6 is a graph of gas velocity crossing a land feature vs. lateral position along 

%1 the land feature for each of a number of apex hal^angles, 9, in accordance with 

™ equation 4, below. 

Fig. 7 depicts a CFD model of the comparative flow field of Example 5C, which 
y is identical to the comparative flow field of Example IC except for the addition of 

J3 20 oxygen consumption, with calculated oxygen flux represented according to the 
indicated gray scale. 

Fig. 8 depicts a CFD model of the "zig-zag" serpentine flow field according to 
the present invention of Example 6, which is identical to the flow field of Example 2 
except for the addition of uniform oxygen consumption, with calculated oxygen flux 
25 represented according to the indicated gray scale, which is the same gray scale as Fig. 7. 

Fig. 9 is a schematic depiction of a "progressive zig-zag" serpentine flow field 
with varying apex half-angle, 9, according to the present invention, depicting 
parameters discussed in the text. 

Fig. 10 depicts a CFD model, which includes oxygen consumption, of the 
30 "progressive zig-zag" serpentine flow field according to the present invention of 
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Example 7, with calculated oxygen flux represented according to the indicated gray 
scale, which is the same gray scale as Fig. 7. 

Fig. 1 1 depicts a CFD model, which includes oxygen consumption, of the 
"progressive zig-zag" serpentine flow field according to the present invention of 
Example 7, with calculated gas velocities represented according to the indicated gray 
scale. 

Fig. 12 is a graph of oxygen flux crossing a land feature in the Y direction vs. 
position in the Y direction for each of Examples 6, 7 and comparative Example 5C, 
taken at an X position approximately 22% of the width of the flow field from the left 
side of the flow field. 

Fig. 13 is a schematic depiction of a progressive serpentine flow field with 
varying segment spacing according to the present invention. 

Fig, 14 is a schematic depiction of a progressive spiral flow field according to 
the present invention. 

Fig. 15 is a schematic depiction of a progressive spiral flow field according to 

the present invention. 

Fig, 16 is a schematic depiction of a progressive spiral flow field according to 
the present invention. 

Fig. 17 is a schematic depiction of a progressive spiral flow field according to 
the present invention. 

Fig. 18 is a schematic depiction of a "progressive zig-zag" serpentine flow field 
according to the present invention having multiple courses. 

Detailed Description 

The present invention provides flow fields for uniform distribution of fluids or 
their active components or properties to and from a target area. The flow field may be 
embodied in a flow field plate or bipolar plate used for distribution of reactants to, and 
removal of products from, opposite sides of a catalyzed membrane in an 
electrochemical cell such as a fuel cell. 

The flow fields according to the present invention are typically embodied in a 
flow field device, which is typically separated from the target area by a fluid transport 



layer. The flow fields according to the present invention provide more uniform access 
of the fluid or its active component or property to the target area, by providing highly 
uniform lateral flux through the fluid transport layer separating the flow field from the 
target area for the transported fluid. 
5 The flow field according to the present invention may be embodied in a flow 

field device such as a bipolar plate (BPP) for an electrochemical cell such as a fuel cell. 
Electrochemical cells include fuel cells, sensors, electrolyzers, and electrochemical 
reactors. Fuel cells utiUze a fuel such as hydrogen and an oxidizing agent such as 
oxygen to produce an electrical current. The two chemical reactants, i.e., the fuel and 
10 the oxidizing agent, separately react at two isolated electrodes containing catalyst. An 
ion exchange element is located between the electrodes to prevent direct chemical 
P reaction of the two reactants and to conduct ions. In the case of a typical hydrogen fuel 

v3 cell, the ion exchange element is an ion conducting membrane (ICM). The ICM 

conducts protons (H"") from the hydrogen electrode to the oxygen electrode. Electrons 
2 15 follow a separate external electrical path, thereby generating an electric current. The 
H combination of an ICM and electrodes is commonly referred to as a "membrane 

electrode assembly," or MEA, The catalyst electrode material may be coated directly 
'fl on the ICM to form a catalyst-coated membrane, or may be coated on the fluid transport 

%J layer discussed below. 

2 20 In conventional fuel cells, MEA's are arranged in a stack separated by rigid, 

electrically-conductive plates which may be known as bipolar plates (BPP). The 
bipolar plate has one or more fluid-conducting channels engraved, milled, molded or 
stamped in the surface(s) facing the MEA(s). The fluid-conducting channels on one 
side of the plate direct fuel to the anode of one MEA while the channels on the other 
25 side direct oxidant to the cathode of the next MEA in the stack. The bipolar plates 
conduct the electrical current generated in each MEA throughout the stack. As used 
herein, "bipolar plate" should be understood to include the end plates of a stack, which 
perform the functions of the bipolar plate on one side only and serve the first and last 
MEA's of the stack. A stack having a single MEA has only two end plates, which are 
30 both encompassed by the term "bipolar plate" as used herein. 
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An additional fluid transport layer is typically situated between the bipolar plate 
and the active catalytic sites of the MEA. In a fuel cell, this may be referred to as the 
diffuser/current collector (DCC), gas diffusion layer, or electrode backing layer. 
(Although the term DCC is used herein, it has become apparent in light of the present 
5 invention that gas transport by processes other than diffusion are important to the 
operation of the DCC.) The DCC is a part of the MEA. Like the bipolar plate, the 
DCC must conduct fluids and electricity to and from a catalyst surface of the MEA. 
Unlike the bipolar plate, typical DCC's are porous throughout and do not function as 
structural members. The DCC typically comprises carbon fiber paper, non- woven roll 

10 goods, or cloth such as ELAT^^ electrode backing material (E-tek, Lie, Natick, MA), 
typically at a thickness of about 0,4 mm. Another typical material is Toray Carbon 
Paper (Toray Industries, Inc., Tokyo, Japan), typically at a thickness of about 02 mm, 
which may additionally be coated with a carbon particle dispersion. 

Uniform distribution of fuel and oxidants over the catalyst electrodes in a fuel 

15 cell should result in more uniform utilization of the catalyst, resulting in better 

performance, stability and durability. Furthermore, this is expected to result in more 
uniform distribution of current density (Amps/cm^), and waste heat generation. This 
should reduce degradation mechanisms that are thermally related, and lead to better 
durability and longer lifetimes. It is desirable then to find methods to distribute reactant 

20 gases most uniformly over the surface of the target area. 

Without wishing to be bound by theory, it is believed that the partial pressures 
of fuel and oxidants at the surface of the catalyst at any given point in an electrode of a 
fuel cell are directly related to the speed of the lateral flux of the gas in the DCC above. 
Gasses are transported to the catalyst surface by a combination of diffusion and 

25 convection. Closest to the catalyst surface it is to be expected that gasses arrive by 

diffusion, subject to Pick's Law. The rate of transport in this diffusion-dominated zone 
is dependent on the gradient in concentration. Farther away from the catalyst surface, 
gas is transported by a combination of diffusion and convection. Since the gradient in 
the diffusion zone depends on the rate at which gasses are replenished by convection in 

30 the zone above, greater and more uniform lateral flux through the DCC should result in 
greater and more uniform flux of reactant gasses to the catalyst surface. 
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Again without wishing to be bound by theory, the importance of convection as 
well as diffusion can be illustrated by considering the Peclet number for mass transport 
to the center of a land area. The Peclet number is one way of comparing the relative 
importance of convective transport to diffusive transport; it is defined as velocity times 
5 distance divided by the diffusivity. For a value of 1, convection and diffusion 
contribute equally. Using values from the Examples considered below, the Peclet 
number for motion in the Y direction from a channel to the middle of a typical land is: 
(0.1 m/s) • (0.001 m) / (L5e-5 m /s) = 6.6; which indicates that convection tends to 
dominate and that both convection and diffusion are important within the DCC. 

10 In addition, it is believed that greater and more uniform lateral flux in the DCC 

improves water management in a fuel cell. It is known that water accumulation at the 
cathode can lead to a reduction in reaction rate by inhibiting oxygen diffusion. This 
condition is referred to as "flooding." It is believed that the present invention can 
eliminate local areas of high flooding by eliminating low flux areas, since a uniform 

15 lateral flux can uniformly carry away water generated at the cathode. 

In contrast to the sort of flow field depicted in U.S. Pat. No. 5,840,438, where 
the active area of the flow field is densely packed with channels, the flow fields 
according to the present invention advantageously include significant land areas. 
Typically 40% or more of the active area of a flow field according to the present 

20 invention is land area; more typically 50% or more, more typically 60% or more, and 
more typically 70% or more. 

The flow field according to the present invention may comprise any number of 
channels. The flow field according to the present invention is typically composed of a 
single channel, but may alternately be composed of a channel composed of multiple 

25 courses effectively parallel to each other, such as depicted in Fig. 18. The active area of 
the flow field may be any suitable size and shape, including rectangular, square, 
polygonal, circular, elliptical and irregular shapes. Typically the active area 
corresponds to the target area, which may be any suitable size and shape, including 
rectangular, square, polygonal, circular, elliptical and irregular shapes. The active area 

30 may be subdivided into separate zones serving separate portions of the target area and 
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treating each portion as a single target area, served by a single channel or a channel 
composed of multiple courses. 

The flow field channels may have any suitable cross-section, including 
rectangular and sloped-side cross-sections. The channels may comprise micro-flow 
5 channels or microstructured features as disclosed in U.S. Patent Applications S.N. 
09/557,712 and 09/430,568, The channels terminate at one end in an inlet, which is 
typically a single opening but may also be multiple openings, and may open into or 
comprise a manifold. Typically, the channels terminate at another end in an outlet, 
which is typically a single opening but may also be multiple openings, and may open 

10 into or comprise a manifold. Where the entire flow of gas into the flow field is to be 
consumed or removed through the fluid transport layer, no outlet is necessary. 

The fluid distributed by use of the flow field according to the present invention 
may be any fluid, including gasses, liquids, supercritical fluids, or combinations thereof. 
The fluid comprises an active component or has an active property, defined above, hi a 

15 fuel cell, the active component may be the oxygen content of air supplied to the cathode 
side of the cell or the hydrogen content of a reformate fuel gas supplied to the anode 
side of the cell. If a pure hydrogen source is used, the active component is the entire 
fluid. In the case of hydrogen or oxygen, the active component is used up as it is 
partially or completely consumed in an electrochemical reaction. The use rate is the 

20 rate at which the active component is withdrawn from the fluid in the transport layer for 
consumption in the target area. In another application, the flow fields according to the 
present invention may be used for uniform distribution of cooling fluids, in a fuel cell 
or some other device. In that case, the active property is the thermal energy content of 
the coolant. In a fuel cell, the air and fuel flows may themselves function as cooling 

25 fluids. The use rate is the rate at which thermal energy is withdrawn from the active 
area into the fluid in the transport layer. In another application, the flow fields 
according to the present invention may be used for uniform distribution of solvents, in 
which case the active property is the solvating capacity of the fluid and the use rate is 
the rate at which solutes are dissolved from the active area into the fluid in the transport 

30 layer. In another application, the flow fields according to the present invention may be 
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used for uniform distribution of air in a planar air bearing, in which case the active 
component is the entire gas and the use rate is the rate at which air escapes the bearing. 

The flow field according to the present invention is embodied in a flow field 
device. The flow field device may be made of any suitable material, but should be 
5 stable to the fluids transported and the conditions of use. Flow field plates or bipolar 
plates according to the present invention for use in fuel cells are typically made of an 
electrically conductive material such as a metal, including titanium and stainless steel or 
conductive carbon materials such as graphite or carbon composites. Alternately, 
materials may be used which are plated, vacuum coated or otherwise coated with 

10 conductive anti-corrosion layers by wet methods, vacuum methods, or any suitable 
method. The channel or channels of the flow field are cut, molded, stamped or 
otherwise formed into the flow field device by any suitable method. The bipolar plates 
of the present invention may be made by any suitable method, including the methods 
described in U.S. Patent No. 5,728,446, to Johnston, et, aL, and pending U.S. Patent 

15 Applications S.N. 09/099,269, 09/557,712 and 09/430,568, all of which are 
incorporated herein in full by reference. 

In one embodiment of the present invention, a flow field is cut, molded, 
stamped or otherwise formed into the material comprising the flmd transport layer. The 
adjacent '^bipolar plate" may then be substantially flat, or may bear additional flow field 

20 channels. 

The flow field according to the present invention is advantageously used in a 
fluid distribution assembly comprising: i) a flow field device embodying the flow field 
according to the present invention and ii) a fluid transport layer disposed between the 
active area of the flow field and a target area. 

25 The fluid transport layer, defined above, may be the DCC of a fuel cell, 

discussed above. The fluid transport layer may be made of any suitable porous or 
permeable material. For use in fuel cells, the material is typically an electrically 
conductive carbon-based material. The fluid transport layer may be any suitable 
thickness. For use in fuel cells, the fluid transport layer is typically less than 1 nun, 

30 more typically less than 500 microns, and can be less than 300 microns. The fluid 
transport layer is typically at least 50 microns in thickness. The fluid transport layer 
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may have any suitable in-plane permeability for gasses, but is typically not greater than 
1 X lO"^ m^, typically between 1 x 10"^ m^ and 1 x 10~13 m2 ^nd more typically 
between 1 x lO'l^ yj^2 ^nd 1 x lO"!^ xn^. 

The fluid distribution assembly according to the present invention provides 
5 more uniform access of the fluid to the target area, by providing highly uniform lateral 
flux through the fluid transport layer separating the flow field from the target area. For 
purposes of the present invention, fluid flux may be determined by analysis of a flow 
field design using computational fluid dynamic (CFD) calculations. CFD is well 
established as an investigative tool for analyzing complex systems of fluid flow, and 
10 therefore any reliable CFD code may be used. Preferably the CFD code "Fluent" (ver. 
5J 5.5, Fluent, Inc., Lebanon, New Hampshire) is used. This is a general purpose, finite 

volume code which uses a form of the SIMPLE method for computing the cell 
4S pressures and fluid velocities, using the assumptions of steady-state laminar flow with 

f% constant viscosity and temperature. Typically, the gas may be considered ideal. Porous 

^ 15 media, such as the fluid transport layer, are modeled using Darcy's Law, for which the 
p. pressure drop is linearly dependent on the viscosity and flow velocity and inversely 

dependent on permeability. (G.K. Batchelor, An Introduction to Fluid Mechanics, p. 
2f 224 (Cambridge University Press, Cambridge, Great Britain, 1967)), ( Darcy, H. 

M (1856). Les Fontaines Publiques de la Ville de Dijon, p. 590, Victor Dalmont, Paris,) 

20 The fluid distribution assembly according to the present invention typically 

provides uniformity to the extent that lateral flux of the active component varies by no 
more than 35% through at least 90% of all overland portions of the fluid transport layer. 
More typically, lateral flux of the active component varies by no more than 30% 
through at least 90% of all overland portions of the fluid transport layer, more typically 
25 by no more than 25%, and more typically by no more than 20%. It is assumed that 
there is a non-zero flow rate of a fluid into the inlet of the flow field, since, of course, 
any flow field will be uniformly still at zero flow. It is assumed that the desired 
uniformity can be achieved for at least one rate of use of the active component of the 
fluid, e,g. one rate of oxygen consumption, as described more fully below and in the 
30 Examples herein. 
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Any suitable flow field design which provides uniform flux through the 
overland portions of the fluid transport layer can be used. Typically, designs having 
non-parallel sequential channel segments are used. When the rate of use of the active 
component of the fluid is greater than zero, typically progressive designs are used, 
5 where the distance between analogous parts of sequential major segments is 
progressively closer toward one end of a flow field channel, or where land areas 
separating the major segments decrease in size progressively with distance from the 
inlet. 

Designs having non-parallel sequential channel segments may include a "zig- 
10 zag" serpentine designs, comprising at least one serpentine channel having non-parallel 
f ^ sequential major segments, such as depicted in Fig, 3 and schematically in Fig. 2. The 

'1; major segments may be curved, but are typically straight line segments. Turning 

iJl segments may be made up of curved segments or one or more straight segments. 

m Alternately, the major segments may meet at a point. 

15 In a progressive design, the distance between analogous parts of sequential 

^. major segments decreases monotonically with distance from the inlet as measured along 

the channel. Likewise, the size of land areas separating sequential major segments 
n decreases monotonically with distance from the inlet as measured along the channel, 

P specifically, the size measured as geometric area or "areal size". Progressive designs 

20 include "progressive zig-zag" designs, comprising a "zig-zag" design where the spacing 
of adjacent channel segments is progressively closer toward one end of a flow field 
channel, such as depicted in Fig. 10 and schematically in Fig. 9. A parallel-segment 
serpentine might also be made progressive, as depicted in Fig. 13. Progressive spiral 
designs are also contemplated, such as depicted in Figs. 14-17. The designs in Figs. 14- 
25 16 require an outlet that is not in plane with the flow field. 

While CFD models, such as used in the Examples below, can provide a more 
complete representation of flow fields according to the present invention, design 
parameters for a given flow field design may also be determined by application of 
geometry and Darcy's Law, as demonstrated following for the case of a zig-zag flow 
30 field according to the present invention. A fuel cell stack may be designed for a 

preferential operating condition of operating current density and mass flow rates. In 
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this case it is expedient to have the flow field optimized for those operating conditions. 
For a given flow field design and size, the pressure drop between the inlet and outlet 
will vary with total naass flow rate. At the preferred operating conditions, the flow rate 
will be fixed so the overall pressure drop will be fixed for the flow field. 
5 Figure 2 defines the geometry and parameters of a single zig-zag flow field 

according to the present invention. The apex half-angle, number and size of loops and 
the DCC permeabilities can be varied to optimize the tradeoff between pressure drop 
and uniformity and magnitude of the gas velocities flowing over the lands. The 
expression developed here gives a first order guideline to selection of those parameters. 
10 The expression can then be adapted to allow treatment of the case in which there is gas 
consumption. 

Figure 2 defines the geometry and parameters of the zig-zag flow field. A single 
loop has a minimum channel spacing on the end of width, w, and a maximum spacing 
on the other end of with, H. The length of the channel leg between the loop ends is L. 
15 Gas flows primarily in the channel, but also flows over the land, through the DCC, 
which has in-plane permeability, Kl- The coordinate x specifies the point at which we 

want to derive an expression for the velocity in the y-direction, U yi^) ^ across the land 
at point X. From Darcy's law, the velocity is related to the gradient in the gas pressure 
as, 



20 



(1) ^ , 

where |i is the gas viscosity and Kj^ the in-plane permeability of the porous DCC 

material in the X-Y plane. 

Since the gas flow in the channel is generally much larger than that flowing 
25 through the DCC across the land, the assumption is made that the pressure at any point 
across the DCC is determined by the flow in the channel, not the DCC. We also 
assume that the pressure decreases linearly down the channel, so that if there are N 

identical loops and the total pressure drop is ^ " ^'"^^^ " P^^^^^t , then the pressure drop 

around any individual loop is , 
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The gradient of pressure over the land at any value of xq, is thus determined by 
the pressure drop along the loop channel between coordinates (xq, -yc/2) and (xq, yo^2), 
divided by the distance, yo- The pressure drop across the land along y at position xq is 
proportional to the total drop around the loop and the fractional distance between 
coordinates (xq, -y(y2) and (xq, yo/2), measured along the channel, or 



(2) 



w + 2'L 



since / ^ in Fig. 2 above. 



Also from Fig, 2, -^^ ^ '^'^ 6 + w j^^pj^^^j^g 3^0 jxi (2) with this gives, from (1), 



Kr^ AP [''^^^ /cose 
' // 7V(m; + 2L) (vv; + 2xotan0) ^ 

10 As angle 6 goes to zero, (3) reduces to the equivalent expression for the case of parallel 
serpentine channels. 

Equation (3) can be used to calculate the velocity as a function of position x, for 
different values of the parameter, 9 . To use (3), the pressure is in units of Pascal, 

viscosity in kg/m-sec, permeability in units of m2, all dimensions in meters, and the 
15 velocity in m/sec. For the flow field of Example 2 below, at ~ IL/min flow rates the 
pressure drop from inlet to outlet is about 5 psig, or 34,500 Pa. For the narrow end of 
the loop we take w = 0.080" = 0.002 m. As the apex half-angle, 9, changes, the number 
of loops that can fit within a defined electrode area, A, changes. Hence, either N can be 
fixed and the area A (or the rectangular length) of the electrode can be allowed to 
20 change, or the number of loops can change with 0 and the area A can be fixed. If the 

area A is fixed, then referring to Fig. 2, it can be seen that for a side length A 1^2^ the 
number of complete loops N that fit within that side length is the lowest integer value 
of N = a1/2 /2(Lsin0 + w). This value for N can be substituted into equation (3) above 
for this case. However, for simplicity for this example, we consider N to be fixed at N 
25 =21 loops and take L - 0.071 m, the width of an electrode area which if square would 
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have A= 50 cm^. Finally using an in-plane permeability of^r=12xl0'12 
equation (3) becomes: 



(4) 



Uy(x) = 6.52-10- 



/cosd 



(0.002 + 2 . tan d) ^ ^^^^^ ^^^^^ ^^^g^ 



5 Fig. 6 shows a plot of equation (4) for values of ^ = 0, 1, 2, 5, 10 and 15 

degrees. The non-uniformity of overland flow speed drops extremely rapidly as the 
angle of the zig-zag increases. Essentially, the variation is gone, except near the apex 
of the loops, already for 6= 5 degrees. This dramatic improvement in gas flow rate 
uniformity with small angles between the loops, as predicted by equation (3), is fully 
10 consistent with the dramatic change in uniformity seen in the CFD calculations shown 
in example 2 above, where the angle is only - L29 degrees ( ~arctan[(H/2-w/2)/L] = 
arctan[(1.8-.2)/71]). hi fact, the magnitude of the flows, e.g. at 6 > 1 degree of Uy < - 

0,2 m/sec, for such a simple model are quite close to the value of 0.13 m/sec seen in 
the CFD model in Example 2, below, which does not make the simplifying assumption 

15 that channel pressure is independent of flow across the land. 

The magnitude of the gas speed through the DCC can be controlled by its 
permeability, K, as shown by equation (3). Gas speed increases with higher DCC 
permeability. The rate at which the flow speed becomes uniform (within a certain %) 
with increasing 0 depends on the flow field parameters w and L, since they enter into 

20 the coefficients and sums of the cosine and tangent terms in (3). The overall geometry 
of the flow field would also seem to be important. E.g. if it were not square, but had an 
aspect ratio of 4, meaning it was 4 times longer than it was wide, how the flow channels 
zig-zagged over the area would be important. If the area was covered by fewer but 
longer loops running nearly parallel to the long direction, the pressure drop per loop 

25 would be larger than if there were more, shorter loops oriented orthogonally. 

As Fig. 8 indicates, consumption of the active component, e.g. oxygen, can 
create additional non-uniformity. A "progressive zig-zag" flow field according to the 
present invention can be designed to partially compensate for the loss of oxygen partial 
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pressure due to consumption by progressively narrowing the apex angle closer to the 
flow field outlet. 

We now consider the case of having oxygen flux uniform down channel, from 
one loop to the other. This will help maintain a uniform local current density, J(X,Y) 
5 from inlet to outlet since J depends on catalyst potential and the flux of oxygen to the 
surface of the catalyst. This oxygen flux close to the catalyst surface in-tum depends on 
the concentration at the surface and the concentration of oxygen in the DCC closer to 
the flow field, as discussed above. 

The approach is to consider the case where the angle ^of the zig-zag decreases 
10 from the inlet to the outlet. As seen in Fig, 6, the average flow velocities for a given 
angle increase as the angle decreases. We consider the case illustrated in Fig. 9, where 
each loop, / = 1 to / = N, has a different value of the channel apex wj, side length L], 

and angle 

From (3) above, for each loop then the uniform velocity is given by : 



15 



(5) 



Since the objective is to determine the conditions for uniform O2 mass flow from loop 
to loop, we can compare the flow at the same point on each loop, e.g. the value when x 
20 = L/ cos(9/) at the widest part of each loop. It is most useful to use an approximation to 

(5) for the following. Since 0 / is small cos(0/) - 1 and Lj**cos(0/) = for all loops. 
Also, at X = L^-, (w/ + L/ ) = L^* . So (5) reduces to: 

(6) U , 

which is true when 2 sin 6 i » wf , as is the case for reasonable values of these three 
25 parameters. This shows that for down-stream loop to loop uniformity, varying G { will 
be the most effective variable to change. 
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We have asserted above that as a first step in keeping the oxygen partial 
pressure over the catalyst surface uniform down channel, keeping the mass flow rate 
uniform through the DCC above the catalyst is useful. The differential mass flow rate 
of oxygen per unit area over the land at the position x = L/ is the flux of oxygen 

entering at that position and is simply the product of the velocity and density of oxygen. 



at that point of loop /: 



(7) ^ {x^L)= ^ " , 

dA RTju 2L.sin^. 

where is the molecular weight of oxygen, R the gas constant, and dA = hdL is the 

element of area through which the gas is flowing in the DCC of thickness h. 

10 Expressions for AP/ and ^'-^^ can be obtained as follows. We assumed above 

equation (2) that the total pressure varied linearly along the flow channel, so the 
pressure drop is simply proportional to the path length. We assume this is true for the 
total pressure when current is generated also. The path length around the i*^ loop is 2L( 
+ w/ ~ 2 Lf. If the total flow field channel length is C, then 

(2L. 4-W-) 2L. 
C C 

For ^'-^^ , we note that the oxygen partial pressure decreases as the loop index 
increases due to consumption. We have to assume some relationship between the two, 
so for simplicity will assume a linear relationship so that if Pin(02) ^he inlet oxygen 
partial pressure, the oxygen is removed in direct proportion to the distance down the 
20 channel, or in direct proportion to the loop index, i. If r| is the stoichiometric flow 
ratio, then l/r\ is the oxygen utilization factor and so 



(9) 



Using (8) and (9) in (7) gives finally. 
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dm,o -K,M Pi„iO^)AP ( i ^ 

(10) '■^{X = L)= ^ ^ 'n^ 1 1_\ 

dA RTfi Csin^,. \ Nr] j 

Uniform oxygen mass flux in all the loops, e.g. particularly the first and last ( i= 
1 and N), can be obtained by requiring the right side of (10) to be the same for both 
values of i, or : 

sin 0, _ (l - 1 A/V^) _ T} 

(11) 



sin^, (l-l/Tj) (^-1), since N»l. 



So in order to have a uniform mass flow of oxygen through the DCC of each 
loop, to a first approximation the angle of the zig-zag channel should decrease as the 
loop index increases towards the gas exit. For the general loop index (/), in terms of the 
initial loop angle Gj, 

si„^,=sin^,|^>. 
10 (12) (^^-1) 



Since r| =2 is a typical desired value for the air stoichiometry of a fuel cell, if the 
loop angle was 5 degrees for the first loop, as Fig. 6 would indicate gives good 
uniformity, then the last loop should have a loop angle smaller by sin 9^ = Vi sin 5*" 
=0.0436, implying 0^ = 2.5"^. Going to a smaller angle would tend to decrease the 

15 uniformity within the last loop compared to the first loop, as indicated by Fig. 6, so 

some compromise between intra-loop and inter-loop uniformity is probably necessary, 
but this example shows that such a trading off is possible to potentially obtain a more 
uniform overall current density than can be obtained with simple serpentine flow fields 
that have parallel and uniformly spaced channels. 

20 This invention is useful in constructing electrochemical cells such as fuel cells. 

Objects and advantages of this invention are further illustrated by the following 
examples, but the particular materials and amounts thereof recited in these examples, as 
well as other conditions and details, should not be construed to unduly limit this 
invention. 



Examples 

In all of the following examples, flow field designs were analyzed using 
"Fluent" (ver. 5.5, Fluent, Lie, Lebanon, New Hampshire) for computational fluid 
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dynamic (CFD) calculations. This is a general purpose, finite volume code which uses 
a form of the SIMPLE method for computing the cell pressures and velocities. 

The gas densities and viscosities used modeled those encountered on the air 
(cathode) side of a fuel cell, where an air flow provides oxygen to the celL Similar 
5 models for the anode side could be constructed using gas densities and viscosities for a 
hydrogen fuel flow. The assumptions used in these models were steady-state laminar 
flow with constant viscosity (2.1 x 10'^ kg/m-s) and temperature (75°C); no reactions or 
gas consumption were included except oxygen consumption as noted. Models that 
include oxygen consumption assume a constant and uniform oxygen consumption rate 
10 over the target area. The gas was considered ideal, implying that the density changes 
directly with the pressure. The porous medium of the fluid transport layer was modeled 
R using Darcy's Law, for which the pressure drop has first order dependence on the 

llJ viscosity, permeability, and flow velocity. In each case, air flow into the flow field inlet 

^ was 800 seem and exit pressure was 253,000 Pascals (2.5 atmospheres). 

I 15 

fi Example IC; Parallel-Segment Serpentine Flow Field. Zero Consumption 

f ^ (Comparative) 

jO A CFD model was built for a comparative flow field, having a single serpentine 

y channel composed of parallel adjacent segments, for a 50 cm^ active area fuel cell, 

P 20 consisting of 25 parallel segments (making 24 loops) each 0.8 mm wide and 1 .0 mm 
deep and spaced 2.0 mm apart (the land width). A CFD model was built in which this 
flow field was in intimate contact with a DCC having a thickness of 122 microns (4.8 
mils) and an in-plane permeability ofK=12xlO-12 ^2 This value of K was chosen 
to represent a typical value for a DCC in an operating fuel cell. The CFD model 
25 calculations were done for the case of 800 seem air flow, and no oxygen consumption. 
Fig. 1 depicts the comparative flow field and shows, in the indicated gray scale, 
calculated gas velocities flowing over the lands in the Y direction within the DCC. 
Velocities are distributed along a spectrum from the Ughtest areas, representing 
velocities greater than 0.2 m/sec, to the dark areas, representing velocities less than 1/4^^ 
30 that of the lightest area. From one side of each loop of the flow field channel to the 
other, there is considerable nonuniformity in the gas velocity over the land area. 
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Examples 2-4; ^^Zig-Zag'' Channel Serpentine Flow Fields. Zero Consumption 

In Examples 2-4, CFD models were built for flow fields according to the present 
invention where adjacent channel segments of a serpentine flow field were slightly 

5 angled with respect to one another, so that any given loop of a channel has an apex half- 
angle, 9, as defined in Fig. 2. To remain suitably comparable to Example 1, the flow 
field of Example 2 also consisted of 25 channel segments each 0.8 mm wide x LO mm 
deep. The active area is 50 cm^ and the overall channel length was approximately the 
same as in Example 1 . Since the 25 channels are 7 1 mm long and 0.8 mm wide, they 

10 cover 14.2 cm^, which is 28.4% of the 50 cm^ active area, leaving 71.6% land area. 

The same DCC was used, having a thickness of 122 microns (4.8 mils) and an in-plane 
permeability of K = 12 x 10-12 ^2 The total air mass flow rate was 800 seem. 

For Example 2, the apex half-angle was approximately L29 degrees and the 
land widths varied from 0.4 mm to 3.6 mm. Fig. 3 depicts the flow field and shows 

15 calculated gas velocities flowing in the Y direction over the lands within the DCC in 
the same gray scale as Fig. 1. For approximately 90% of the overland active area, the 
velocities are approximately 0.13 m/sec. Some higher velocities are seen near the inlet 
and outlet holes. This represents a large improvement in uniformity over Example IC. 
Fig. 4 is a graph of inlet-to-outlet pressure drop as a function of mass flow rate 

20 for the Example 2 and comparative Example 1 flow fields. It is seen that they have the 
same pressure drop, and therefore pressure drop was not disrupted by the change from a 
parallel-segment serpentine to a zig-zag serpentine design. 

Examples 3 and 4 were similar to Example 2 except that apex half-angles 9 
were changed to lower and higher values, 0.89 and 2,10 degrees respectively. The 

25 separation between sequential segments of the serpentine channels at the widest part 
(H) was 2.6 mm, 3.6 mm and 5.6 mm respectively for Examples 3, 2 and 4, Due to the 
change in apex half-angle, the flow field of Example 3 had 31 segments (30 loops) and 
the flow field of Example 4 had 19 segments (18 loops), and therefore the total channel 
length and pressure drop were larger for Example 3 and less for Example 4. It also 

30 follows that land area was less for Example 3 and greater for Example 4. Fig. 5 is a 
graph of gas velocity in the Y direction across a typical land feature as a function of 
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lateral position along the land feature (i.e., in the X direction as depicted in Figs. 1 and 
3), for each of Examples 2-4 and comparative Example 1. Gas velocities in each of 
Examples 2-4 are more uniform than in comparative Example 1, although Example 2 
appears to be more optimal than Examples 3 and 4. 
5 The plots in Fig. 5 also show that the magnitude of the average velocities is 

reduced as the apex half-angle increases. As shown in more detail above, this effect 
can be offset by changing the value of the permeability, K, of the DCC. 

Examples 5C and 6: Non-Zero Consumption 

10 In Example 5C, the CFD model for the comparative flow field of Example IC, 

having a single serpentine channel composed of parallel adjacent channels, was run 
with the modification that oxygen consumption was added to the model at a 
stoichiometric flow ratio, r|, of 2.0. The stoichiometric flow ratio is the ratio of oxygen 
supplied to oxygen consumed. Fig, 7 depicts the resulting in-plane oxygen flux, which 

15 is predominantly in the Y direction. 

In Example 6, the CFD model for the flow field of Example 2, having a single 
"zig-zag" serpentine channel, was run with the modification that oxygen consumption 
was added to the model at a stoichiometric flow ratio, r|, of 2,0. Fig, 8 depicts the 
resulting in-plane oxygen flux, which is predominantly in the Y direction, 

20 Figs. 7 and 8 demonstrate that non-unifomaities are amplified when oxygen 

consumption is initiated. The "zig-zag" serpentine is far superior in uniformity. 

Example 7; ' Vrogressive Zig-Zag'" Serpentine Flow Field 

A CFD model was built for a flow field similar to Example 6, including oxygen 
25 consumption at a stoichiometric flow ratio (r|) of 2.0, but differing in that the apex half- 
angle, 6, as defined in Fig. 9, is progressively narrowed from 1.75 degrees at the inlet to 
0.80 degrees at the outlet. The small end of the lands is fixed at 0.4 mm while the large 
end varies from 4.8 mm to 2.4 mm. Note that the average large end is 3.6 mm, 
matching Example 2, and the number of channel segments is 25 as well, also matching 
30 Example 2. This "progressive zig-zig" flow field also comprised channel segments 
each 0.8 mm wide x 1.0 mm deep, forming an active area of 50 cm^ and the overall 
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channel length and land area were approximately the same as in Example 2. The same 

DCC having a thickness of 122 microns (4.8 mils) and an in-plane permeability of K = 

12 X 10-12 jn2^ and 800 seem total air mass flow rates was used. 

Figs. 10 and 1 1 demonstrate the in-plane oxygen flux in the Y direction and the 
5 gas velocities in the Y direction for this progressive zig-zag flow field. The 

progressively narrowing apex angle causes overland gas velocity on the downstream 

side of the flow field to increase, thereby compensating for the progressively decreasing 

partial pressure of oxygen caused by consumption. 

Fig. 12 compares the flow fields of Examples 5C, 6 and 7. It is evident that the 
10 zig-zag serpentine flow field according to the present invention demonstrates far more 

uniform oxygen flux than the comparative parallel-segment serpentine, and the 

progressive zig-zag demonstrates even better uniformity. 

Various modifications and alterations of this invention will become apparent to 
15 those skilled in the art without departing from the scope and principles of this 

invention, and it should be understood that this invention is not to be unduly limited to 
the illustrative embodiments set forth hereinabove. All publications and patents are 
herein incorporated by reference to the same extent as if each individual publication or 
patent was specifically and individually indicated to be incorporated by reference. 

20 
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